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INTRODUCTION 

Lobe mixers are devices which are currently being used for thrust augmenta- 
tion on a variety of turbofan engines. The mixer is designed to provide a means 
of generating mixing between the turbine and fan streams before they enter the 
exhaust nozzle. Through careful design, the mixing can result in more uniform 
energy distribution at the nozzle exit. Provided that the losses incurred in the 
mixing process are minimized relative to the amount of mixing, small but signxfi 

cant performance gains may be realized. 

Obviously, due to the sensitivity of mixer performance to changes in lobe 
shape and number, inlet velocity and temperature ratios, and secondary flows 
generated by the lobes, a computation procedure which could be used to sift 
through a matrix of possible design configurations would be of significant value. 
Such a procedure was previously developed by Krestovsky, Briley and McDonald 
(Ref. 1), based on the velocity decomposition approach of Briley and McDonald 
(Ref. 2). The analysis synthesizes concepts from potential flow theory, 
secondary flow theory, and from an extension of three-dimensional boundary 
layer theory in a manner which allows efficient numerical solution by forward 
marching in space. The resulting procedure was used extensively in Refs. 3-5 
to examine the effects of lobe shape, lobe generated secondary flows, and 
turbulence effects on mixing. The results obtained in Refs. 3-5 were, in 
general, found to be in excellent agreement with available experimental data. 

Through interpretation of the results obtained using the analysis of 
Ref. 1, a more detailed understanding of the mixing processes in lobe mixers 
was obtained. As a result, it became apparent that it would be beneficial to 
extend the analysis of Ref. 1 and further develop and improve its capabilities. 
These extended capabilities, which are the subject of the present investigation, 
allow the current version of the analysis to be applied to mixer flows with 
swirl, and to allow more general specification of initial condition with 
nonzero secondary flow vorticity. In addition, the procedure for determining 
the required a priori compressible potential flow pressure gradients from an 
incompressible flow solution, has been improved. Finally, the formulation of 
the k— e turbulence model used in the previous version of the procedure was 
somewhat unreliable, leading in some instances to predictions of negative 
values of turbulence kinetic energy. Under the current investigation, the 
reliability of this model was improved by making minor changes in the numerical 
formulation, and optional turbulence models were added to the procedure. 
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A. a result the present version of the dr analysis end conputer code 
represents a reliable, proved analysis uith extended capabilities uhich can he 
:i CO further the understand of lohe dr flow fields, and to aid in 

mixer design. 



ANALYSIS 
General Approach 


The present analysis is described in detail in Ref. 1, and is outlined 
only briefly here. The analysis is based on approximations made relative to 
a curvilinear but not necessarily orthogonal coordinate system fitted to and 
aligned with the flow geometry under consideration (cf. Fig. 1). The coordinate 
system is chosen such that the streamwise or marching coordinate downstream of 
the lobes either coincides with or is at least approximately aligned with a 
known inviscid primary flow direction, as determined for example by a potential 
flow for the given geometry. Transverse coordinate surfaces must be either 
perpendicular or nearly perpendicular to solid walls or bounding surfaces, 
since diffusion is permitted only in these transverse coordinate surfaces. 

Equations governing streamwise vorticity and a scalar viscous correction 
o to a known inviscid primary flow velocity \ are derived utilizing assump- 
tions which permit forward-marching solution, provided reversal of the composite 
streamwise velocity does not occur. Terms representing diffusion normal to 
transverse coordinate surfaces are neglected. Approximate pressure gradients 
are derived from the inviscid primary flow and imposed in the streamwise momentum 
equation. Secondary flow velocities are determined by scalar and vector poten- 
tials in the transverse coordinate surfaces, to suppress streaimrise elliptic 
behavior requiring downstream boundary conditions. Finally, the use of stream- 
wise vorticity in obtaining secondary flow velocities avoids the explicit use of 
transverse momentum equations in predicting the secondary flow velocities. 

Primary- Secondary Velocity Decomposition 

The analysis is based on decomposition of the overall velocity vector 
field U into a primary flow velocity U p and a secondary flow velocity IK The 
overall or composite velocity is determined from the super-position 


For lobe mixer flows the primary flow velocity is represented as 

( 2 . 2 ) 

Up “ u i u V 

where Uj is a known inviscid primary flow velocity satisfying slip conditions 
and determined from an a priori potential flow solution in the geometry' under 
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consideration. The (non-dimensional) scalar quantity u y is a viscous velocity 
profile factor which introduces viscous shear layers and may also correct for 
internal flow blockage effects. In the present case of the lobe mixer u y accounts 
for both boundary layers and the free shear layer which exists between the fan 
and turbine streams. The viscous velocity correction u v is determined from 
solution of a primary flow momentum equation. As will become apparent, for the 
lobe mixer it is conventient to solve for the component of U p in the primary 
flow direction rather than u^. 

To suppress streamwise elliptic behavior, the secondary flow velocity U, 
is defined by, and presumed derivable from, scalar and vector surface potentials 
denoted <j> and i|», respectively. If ij denotes the unit vector normal to trans- 
verse coordinate surfaces (also presumed here to be in the direction of the 
marching coordinate), if P is density, and if P q is an arbitrary constant 
reference density, then Ug is defined by 

u s = +(/>„//>) Vx (2 ' 3) 


where V is the surface gradient operator defined by 
s 

v, H V - T t ( i • V) 


(2.4) 


It follows that i^ • U s = ° and thus Ug lies entirely within transverse 

coordinate surfaces. Equation (2.3) is a general form permitting both rotational 

and irrotational secondary flows and will lead to governing equations which are 

ip tic in transverse coordinate surfaces and which are, therefore, solvable 

within a forward-marching context. The overall velocity decomposition (2.1) 

can be written _ _ — , # , * « » 

U - U T U +V <p + (/> 0 / />) Vx, ,r 

1 * (2.5) 

-U D +V # * + (/>.//») vxty 


Surface Potential Equations 

Equations relating <{> and with u^, p, and the streamwise vorticity 
component ^ can be derived from Eq. (2.5) as follows: 

V~/>U - O - V-/>U,u v +V-pV s <f> +/> 0 V-VX? ( ^ 
l*VXU s XI • l.*VXU I U v + i • Vx(a //>)VX I + i ( - VXV # ^> 


( 2 . 6 ) 

(2.7) 
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Sin<t- the last term in each of Eqs. (2.6, 2.7) is zero by vector identity, 

Kqs. (2.6, 2.7) can be written as 

V-pVg4> " ( 2 . 8 ) 

V vx( / > o / / ,,vx ty " fl « -»r vx °i u v (2 * 9) 

which are elliptic for $ and tp in transverse coordinate surfaces. The last 
term in Eq. (2.9) is identically zero in a potential flow coordinate system 
for which ij and Uj have the same direction, and is small if ± 1 and U I are 
approximately aligned. Given a knowledge of u v , P., and p, the surface potentials 
4 and ip can be determined by a two-dimensional elliptic calculation in trans- 
verse coordinate surfaces at each strearawise location. In turn, U g can be 
computed from Eq. (2.3), and the composite velocity 0 will satisfy continuity. 
Equations for u v and ^ are obtained from the equations governing momentum 
and vorticity, respectively. 

Primary Momentum and Pressure Approximation 

The streamwise momentum equation is given by 

A r i * _ (2.10) 

^•[(U-VJU + (Vp)//>J - lg - F 

where p is pressure and pF is force due to viscous stress. Terms in F represent- 
ing streamwise diffusion are neglected; however, since the viscous terms are 
complex for compressible flow, the modified viscous force is temporarily denoted 
F’, and further consideration of viscous terms is deferred to a later section. 

The remaining assumption for Eq. (2.10) concerns the pressure gradient 
term and is designed to permit numerical solution as an initial value problem. 

An obvious pressure approximation for curved flow geometries is to impose 
pressure gradients from an inviscid potential flow (Briley [6]), and for 
internal flows, to correct these with a mean pressure gradient term depending 
only on the coordinate (Pantankar & Spalding [7]). This approximation is 
appropriate both for flows consisting of an irrotational core region with 
thin shear layers on bounding surfaces, and also for some fully viscous flows. 

If pressure gradients are derived from the inviscid potential velocity 
0 and if a mean viscous pressure correction P v (Xj) is introduced, the pres- 
sure approximation can be written as 
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i, • Vp * i | -[Vp v (x | }-/)7(0 I -D I )/2] 
£ l 1 Vlp I + p v ) 


( 2 . 11 ) 


where p I is the imposed pressure. Typically, for internal flows, p v is 
determined to ensure that an integral mass flux condition is satisfied, such 
as 



constant 


( 2 . 12 ) 


Combining Eqs. (2.10) and (2.11) and setting F = F* provides an equation 
nominally governing u^: 

j | -[(U-V)U +(Vp ¥ )//> -VfUj-U j)/z] (2 13) 

where U is to be written as 

U - I.UjUy + ^(VjUy +v g ) + i 3 (w 1 u v + w,) (2.14) 

where u^, v^, are components of Uj, and v g , w^ are components of U s - 

Secondary Vorticity 

An equation governing in compressible flow can be obtained from an 
approximate application of secondary flow theory. The equation governing the 
growth of vorticity along a streamline for compressible flow with constant 
viscosity is given by Lakshminarayana & Horlock [8] and may be written for 
variable viscosity as 





I 


s-Vx Vp + 

p p* 


i . 

2 5 


(2.15) 


where q is velocity magnitude, s is the unit vector along a streamline, s is 
distance along a streamline, to is vorticity in the direction of the unit 
principal normal vector n of the streamline, for which R is the principal 
radius of curvature. These quantities are related by the Frenet formula 
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n/R = 3s/Ds = s • Vs. In Eq. (2.15), G " V x F, mid the term contain I ng p 
vanishes if p is constant. Since an intrinsic coordinate system formulation 
as in Eq. (2.15) provides an "inverse" coordinate system if used to compute 
U, and since intrinsic coordinates are degenerate on no-slip surfaces and 
nonor thogonal for general rotational flows, intrinsic coordinates are not 
attractive for numerical computation. However, if the coordinate system used 
for computation is approximately aligned with the flow direction, then an 
approximate equation governing ft, can be derived from Eq. (2.15) by replacing 
s by i, ns in the following development: 


n s s-fi ij-xi xij 

q * s^u i~ U Uj 


(2.16a) 


n n (s.Vs).fl (T^vT,) n Xl ni (2.16b) 

qR s-u ^ “ U( R, 


where u^ = u^t^; is the principal radius of curvature of the coordinate, 
>™d ^nl is vorticit y in tlie direction of n^, the principal normal of the Xj 
coordinate line. The quantities ^ and are defined by the Frenet formula 


n l /R l = i l 


Vi^. To illustrate, in an orthogonal coordinate system. 


A 


n 

i 






+ A_£M 

h.h 3 ax 3 J 


(2.17) 


where h^, h 2 , h^ denote metric coefficients. If p varies, p is replaced by 
the imposed pressure Pj as defined in Eq. (2.11). Finally, taking q 2 = Uj 2 
and s = i^ in the last two terms in Eq. (2.15), and neglecting stre;imrise 
diffusion, Eq. (2.15) becomes 


U-V 



2H 


Bi 





A 





(2.18) 


where G' does not contain streamwise diffusion. 

The transverse vorticity ft^. in Eq. (2.18) contains components which, in 
ortliogonal coordinates and assuming an irrotational U^, arc given by 
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Cl t - ? t *Vxu 




dx, 


h,h 3 dx, 


3r ,h > w *’ 


(2.19a) 




• V 7x0 


^ + it .m* 

h 2 ax t h, ax, 


1 d 

+ mTK' < "* , « > 


it 


(2.19b) 


In the applications contemplated here, the first term on the right-hand side of 
each of Eqs. (2.19a-b) is expected to dominate, and the remaining terms may be 
neglected as a convenience. 

Energy Equation 

Since the flow being considered is nonadiabatic, solution of an equation 
governing energy is required. The energy equation can be written in a variety 
of forms, one of vhich is 


/jU-VE • V-k7T + U-F + 4> (2.20) 

where E is total enthalpy, T is temperature, k is thermal conductivity, and 
* is the dissipation function. Solution of Eq. (2.20) by forward marching 
integration requires only that terms representing streamwise conduction of heat 
and also streamwise viscous diffusion in F be neglected. 

Compressibility Relations 

In the imposition of streamwise pressure gradients, compressibility effects 
are represented by introducing the perfect gas equation of state p = pRT. 

No other assumptions are necessary for consideration of compressible flow. 

For moderate subsonic Mach numbers, inviscid pressure gradients would ideally 
be obtained from a compressible potential flow calculation or otherwise from an 
incompressible potential flow corrected for compressibility using a suitable 
compressibility correction. Replacing p in Eq. (2.11) by the state equation and 
eliminating temperature using the temperature-enthalpy relation 

U-0 (2.21) 


8 



ORIGINAL PAGE IS 
OF POOR QUALITY 


where «: denotes specific heat, the following auxiliary equation relating 
tlu imposed pressure gradients with density, velocity, and total enthalpy 
is ontai**£d: 


V v [ p i + p v (x » } ] ' 



yH 

r 



( 2 . 22 ) 


where y is specific tieat ratio. A slight simpl , fication rcsulLs if U * U is 
— 2 

replaced by (i^, • U) Eq. (2.22). 

Governing System of Equations 

A complete system of six coupled equations governing u^, 5^, 4s <S E, 
and i> is given by Eqs. (2.8), (2.9), (2.J3), (2.20), (2.22), and (2.18). 
Ancillary relations are given . ' Eq. (2.5) for composite velocity, Eq. (2.12) 
for mass flux, and E<3 s* (2.19a— b) for transverse vorticity. 

The Turbulence Models 

The present analysis is to be used to predict a turbulent flow field, 
thus the stress terms F* and G 1 in Eqs. (2.13), (2.18), and (2.20) contain 
turbulent shear terms in the form of fluctuating velocity and velocity- 
temperature correlations. These terms may be modeled using any one of a number 
of turbulence models which have been included in the present analysis. In 
general these turbulent shear terms are modeled through the introduction of a 
turbulent eddy viscosity, and this eddy viscosity may be computed from one of 
three optional turbulence models. The most complex and general turbulence 
model used in the present analysis is the k-c turbulence model described 
by launder and Spalding (Ref. 9). 

The transport equations governing k and c are given in vector form 
(Ref. 9) for steady flow as 


iD U-Vk= V(/i+/i T /o- k )Vk + P- pc 

e 

pU.V«;V-(/x+/i T /«7- € )Vk + C,— P-C 2 

r k k 

In Eqs. (2.23, 2.24), P is the turbulence production. The quantities o^, 

O , C , C are empirical constants which, as recommended by Uiunder & Spalding 
(Ref. 9), are taken as 1.0, 1.3, 1.44, and 1.92, respectively. Upon solui io:, 


(2.23) 


(2.24) 


iT 
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of Eqs. (2.23 and 2.24) the turbulent viscosity and length scale may be 
computed as 

fij = C^pk 2 /* (2.25) 

and 

C/1 w4 k 3/2 A (2 - 26) 

If low Reynolds number flows are considered (i.e. wall boundary layers are 
resolved) the constants and are adjusted from their high Reynolds number 
values. As suggested by Launder and Spalding (Ref. 9), C 2 is given as 

C 2 = C 2 ^(l.0-0.3e" R t) (2.27) 

where is the high Reynolds number value of given above, and is a 
turbulence Reynolds number 

R t = (2.28) 


The constant is adjusted following Shararoth and Gibeling (Ref. Id) 

= 40, 2 


(2.29) 


where a^ is a function of a different turbulent Reynolds number R^, and is 
given by McDonald and Fish (Ref. 11) as 


°,-°o[ i w]/| l - 0 + G - 6G °o[-T5T-']| 


In Eq. (2.30) a is taken as 0.0115 and f(R ) is given as (Ref. 11) 
o x 


f(R r > ■ 100 R r 


0.22 


r t <i 


(2.30) 


(2.31a) 


f(R r ) =66 IR t + 614.3 R T >40 


(2.31b) 


Between the limits 1 < R^ < 40 a cubic is used in the two functional forms 


given by Eq. (2.31). For these functions R^ is defined as 

M I X 




(2.32) 


Although the k-e turbulence model given above provides a general 
turbulence model in which both the velocity and length scales of turbulence 
are predicted, the formulation has proven to be unreliable and unduly 
sensitive to the initial profiles at times leading to the prediction of 
negative values of k and/or c. To provide a means of obtaining meaningful 
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results when such a failure occurs an alternative turbulence model may be used. 
In this turbulence model a free stream length scale is specified and the 
turbulence dissipation equation <Eq. 2.24) is not solved. Instead, the 
turbulence dissipation appearing in (Eq. 2.23) is eliminated using (Eq. 2.26) 
in conjunction with the specified length scale. In near wall regions the 
length scale is damped following McDonald and Fish (Ref. 11) 

1 = i^lonM k 6/J.^ (2 ' 33 

where £ is the freestream length scale, * is the von Karman constant 

co 

(=0.43), and ^ is a sublayer damping function given as 


J3 -- P ,/2 [ (d* - 23) /8] 


Here, P is the normal probability function and d + - d^/v where u t is the 
friction velocity. 

For situations where it might be desirable to investigate the role of 
turbulence models, a third and very simple turbulence model has been included 
in the analysis. This turbulence model is referred to as a wake turbulence 
m( ,dcl and is derived from the definition of turbulent viscosity and 

Prandtl’s suggestion that the Jiffusional flux of some property 

may be expressed - 

... _ 1 it 1,2 f — — = -y-L- (2 35 

" v $ ' h L dy P 1 

hence one obtains the Prandtl-Kolmogorov relationship 

u = — k" 2 L p (2-36 

The wake turbulence model is obtained by assuming that the turbulent length 
scale and velocity scales are constant throughout the flow field. From^ 
equations 2.25 and 2.26 it follows that the constant l/o^ is simply <y 
The velocity scale is taken as the square root of the initially specified 
turbulence kinetic energy k. The length scale is specified using the 
geometric constraints of the mixer as a guide, for instance assuming that the 
initial length scale would not be smaller than the inlet boundary layer or 
larger than the duct height. As a consequence of this formulation, the 
turbulent viscosity is dependent only upon the initial choices of k and £, 
and the local density which is determined as the solution progresses. 
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application to flow in a lobe mixer 


In this section, specific details arising in the application of the 
foregoing analysis to lobe mixer problems are given. A cross section of a 
typical lobe mixer geometry is shown in Fig. 2. The area immediately down- 
stream of the nozzle plug is faired in with an assumed streamline to model 
the reversed flow region expected in this region. Since the flow area thus 
excluded from consideration is very small, this treatment is not believed 
to introduce significant error. A more detailed analysis of this separated 
flow region could be performed following a zone embedding approach using the 
Navier-Stokes equations (cf. [12]). Although the mixer geometry is axisym- 
metric, the flow is three-dimensional due to the aximuthal variation of the 
hot and cold streams. However, due to observed symmetry, only a pie-shaped 
region of the transverse coordinate surface need be considered. The shape 
of this region and the extent of typical hot and cold streams at the mixer 

exit surface are shown in Fig. 3. 

Coordinate System 

Curvilinear orthogonal coordinates x, y, z are constructed to fit the 
flow passage boundaries downstream of the lobes as shown in Fig. 2. Metric 
coefficients h., h , h 3 are defined such that incremental distance s is 
determined by («x>* = (h^x) 2 + (h 2 dy) 2 + ^Sz) 2 . In planes of constant z, 
(azimuth), orthogonal streamline and velocity potential lines from a two- 
dimensional planar incompressible potential flow analysis are utilized as the 
coordinate lines for constant y and x, respectively. This x-y coordinate 
system is then rotated about an axis representing the centerline of the lobe 
mixer, to form an axisymmetric coordinate system. Thus, the z direction can 
be regarded as cylindrical (i.e., 0) and h 3 = r. The x coordinate is taken 
as the primary flow or axial coordinate, and is associated with surfaces for 
which the two-dimensional potential is constant. The y and z coordinates 
define transverse secondary flow planes at any given x location. In this^ 
coordinate system, the normals \ to the transverse planes coincide with ± 1 
the unit vector in the x direction. A two-dimensional incompressible 
potential flow analysis and computer program developed by Anderson [12] was 
employed in the present investigation, without modification, to compute the 
necessary coordinate data. 
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The Imposed Pressure Field 

. . Wnf i the axisymnetric coordinate directions generated 

AS L g C p ^ 1 two-d imens iona 1 potential flow coordinates about an axis 
by rotating pi- • ^ ry ^ secondary flow directions, 

provides 3 ^ r °^. roensional potential flow pressure gradients are 

i; - , r 

»n for compressibility effects ns 

tbo given coordinate •> compressibility correct ion 

- - - r ar4 ~ « 

'** l'" rr ""' C ' i ted i„ the correction. This inaccuracy in 

„ ic h subsonic Mach number „ ln int ercst in 

, „,»«,,re field was accepted xn K_t. 

‘" C T ,Ld H te resl o„ downstream of the iobes and upstream of the 

mixing was the tmmed . b bl „ „ a s not a major issue. Under 

throat, and consequently scali „„ based on the 

‘"VTnlhir:: developed. Civen the incompressible 

work of I.i obi c- (Ref 15) demonstrated that an 

, . ,i i< licblein and Stockman V^ex. x-« > 

lllacl'af the compressible — neld — - obtained 

from the relationship (V°l) 


u c = u i(t-) 


( 3 . 1 ) 


, ro „ is the compressible velocity, p, is the incompressible (or stagnat.on, 
s t / is «- local compressible density, and U, Is the atea averaged 

^ romprcssibl e velocity .at a given axial location deter.rned as 

J u I dA 


int 


u i = 


( 3 . 2 ) 


JdA 


rhat the density ratio, P T /P. be computed using 
Licblein and relat ionshi P s with the Mach number determined 

one-d troens tonal lsen trap appll „ tio „ t h. density rntio is compnt.d 

by the ore. r :"; nships ln conjunction with the local Mach number 

by using d -;;;”7;" foruacd Mrchlpg in this M nncr, a fully 

computed as part velo city field can be approximated from the 

three-dimensional eompres ^ flpp K „„„l.dge of the compressible 

two-dimensional incompressi P .dnressureinthcpolenti.nl 

vdocity, and assuming constant total temperature and pressure 
flow the static temperature is obtained from 
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T 0= T — l/2Ll| 


(3.3) 


The corrected pressure is then found from the relationship 

r 

P=P 0 (T/T c ) r -' < 3 - 4 > 

from which the compressible imposed pressure gradients are determined. 

Governing Equations in Orthogonal Coordinates 

As in Ref. 1, the governing equations reviewed in the previous section, 
and the turbulence model equations must be expressed in the orthogonal 
coordinate system used to represent the mixer passage geometry. Although these 
equations are fundamentally the same as those presented in Ref. 1, the approxi- 
mations to the stress terms in the primary momentum equation, and the production 
term in the turbulence kinetic energy and dissipation rate equations have been 
altered in an effort to improve the reliability of the k-e turbulence model. 

With the exception of these approximations, the equation in orthogonal coordinates 
follows directly from Eqs. (2.8, 2.9, 2.13, 2.18, 2.20, 2.22, 2.23 and 2.24) by 
performing the vector operations in the given coordinate system. 

Throughout the remaining discussion, all variables in the governing equa- 
tions are nondimensional, having been normalized by the ^ - 'wing reference 
quantities: distance, L r ; velocity, U r ; density, P r ; temperature, 7 r ; total 

enthalpy, U^; pressure, p^lT; viscosity, u r - Here the subscript r denotes a 
reference quantity. This normalization lends to the following nondimensional 
parameters: Mach number, M; Reynolds number. Re; Prandtl number, Pr; and 

specific heat ratio, y. These parameters are defined by 

M * U f /c , Re - /> r U f L r //t r , Pr - c p /i r /k , y - C p / C y (3.5) 

where is the molecular viscosity, k :* ; thermal conductivity, and c p and 
are the specific heats at constant pressure and volume. The reference speed 
of sound, c, is defined by c 2 - yR g T r , where R g is the gas constant. Since 
the flow is assumed to be turbulent, the dependent variables are taken to be 
the time-averaged quantities in the usual sense; however, third order correla- 
tion as well as second order correlations including the fluctuating density 
have been neglected. The decomposition of the velocity vector thus can be 
written 
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U • I ,u + 7 f v + IjW 


Where and are the unit vectors in the x, y, and z coordinate 

directions, respectively. Since i^ is closely aligned with the primary flow 
direction, Vj is small and is neglected as a minor convenience. This approxi- 
mation is of little consequence since the neglected quantity v^ effectively 
reappears ns part of v s and is thus determined from solution of the continuity 
equation during the forward-marching solution process, instead of being imposed 
a priori from the potential flow. Since Wj is zero by definition, the velocity 
decomposition can be written as 

“■', U p +T £ V . +T 3 W . (3.7) 

Since u^ is the only streamwise component of U the subscript "p" is 
dropped. Under the stated assumptions the streamwise momentum equation can be 
written as 


, du du du dh 9 dh 

pUh ^~d 7 + ^ )h * h 3 77 + /> (W S )h . h 2 77 + />( y s) Uh 3^ - 


dpy(x) 




- nfui l 2 h atl3 + h h df>1 -i- K K dPv(X) 1 d T »>■ h, 

^ s hz sr + 

d r , , du I 2/i. dh [ dw dh,l 2tih dh, 

Tz ^ + ^ y dTj - dF + v dTj - ^ 

r hg / dh 3\ 2 /_^2\ h 3 ~ | 

L h 3 V d x ' ' dx ' h 2 J U 


h, hj d(u/h, ) 
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with the Prandtl number Pr and turbulent Prandtl number Pr T assumed to be 

To obtain the rotational secondary flow components of the velocity, the 
streamvise vorticity equation is written as 


u / ac 

dpu \ 

. v « h . ( PU ±L _ 


w s 


ae 

pU “T 

Tl^ar 

ax J 

O h 2 v f ay 

£ ay / 

0 

h 3 V 

^ az 
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i au h, 

i r ap 

ap 

ap 

. 41 ] 

* - 

"S7 a * 

h,h s 

M a * 

az 

ay 

az J 


Re 


h, a r h 3 dh,(/t+/i.T)C] A. 


where q 2 = u 2 + v 2 + w 2 and the vorticity is related to velocity by 

i t ah a w > ah e v, \ 

{ “ "h"h^ l ay ' az / 


( 3 . 9 ) 


. 0 . 


( 3 . 10 ) 


( 3 . 11 ) 
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The vector potential ^ is determined from 


1 3 j 

\( h 3 \ * W 1 

1 d 

J- 

Sri 

i «MM1 

h 4 h, ~w\ 

p\\} ay J 

dz 

1 dz j 


(3.12) 


and the rotational secondary flow velocities are determined as 

I 

V > ' />h,h 3 dz 

I dth,^) 

' /> h i h e d y 

The scalar potential is governed by 


r /> h , h s 1 d 


d(h c h 5 />u) 

[ dy J dz 

h 3 dz J 

dx 


The velocity components v and w, are found to be 

9 9 

i d 4>s 
i a<#> 5 

~h~~dF 

The secondary flow velocities are thus 

.- v * + v * 

1 , * + w ^, 


(3.13a) 


(3.13b) 


(3.14) 


(3.15a) 

(3.15b) 


(3.16a) 

(3.16b) 


The turbulence model equations must also be expressed in orthogonal coordinates. 
It is assumed that the turbulence equations are tensor invariant [9], and thus 
may be expressed directly in orthogonal coordinates by performing vector opera- 
tions in such a coordinate system. With streamwise diffusion neglected to allow 
forward-marching solution, the result for the turbulence kinetic energy equation 
in steady flow is 
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(3.17) 


and similarly the turbulence dissipation equation is given by 


fUhjhj-^r + P v * h | h 3 + ^ W » h l h *dX 




(3.18) 
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The production term is approximated as 


P = 


it 

Re 


■£<!«] 


(3.19) 


The turbulent viscosity is then 


H-i “ 


c upk 


R« 


and the length scale may be determined as 

3/4 3/2 


t - 


Cp It 


(3.20) 


( 3 . 21 ) 


When the k-1 turbulence model Is used, E„. (3.19) and (3.20) are cmdrined 
to eliminate the appearance of the turbulent dissipation from Eq. (3.19). 

t = c/ /4 >? n /X W ' 22) 
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and it follows that 




C/ 4 Vy-JL Re 


(3.23) 


Specification of Initial Conditions 

The initial conditions for a lobe mixer calculation may be specified by 
using either an automated starting routine or by reading in experimental data. 

The fully automated procedure is considered first. 

To obtain initial conditions for a lobe mixer calculation, it is necessary 
to specify the velocity and total temperature of the respective hot and cold 
streams in addition to a mean value of static pressure. The lobe shape is 
specified, and based on the specified lobe shape, a decision is made as to. 
whether a grid point lies in the hot or cold stream. The appropriate values 
for velocity, temperature and density are then assigned to the grid point. 

The total energy is assumed to be constant but may differ in both streams. 

The velocity in the two respective streams is assigned its nondimensional 
reference value which is then corrected to account for normal pressure 
gradients present at the initial plane as determined from the axisymmetric 
potential flow. To account for boundary layers on the lobe, hub and shroud 
surfaces, the free stream velocity profiles are further scaled in accordance 
with an assumed turbulent boundary layer velocity profile, and the distance 
from the surface. Radial secondary flow velocities at the initial station 
may be specified parametrically in each stream as a fraction of the streamwise 
velocity. This in effect sets the flow angle in the fan and turbine streams, 
respectively. Finally, a swirl velocity component may be specified, when 
appropriate. 

The initial velocity profile with boundary layers is given by 


u(u 


STHM 


» s «»y /S BL» s * hJ 


(3.24) 


where U is the velocity of the appropriate stream, 6 gL is the boundary 
layer thinness, y/6 represents the shortest distance from a surface (hub, 
shroud or lobe) to the grid point in question, 6 is the local displacement 
thickness and H is the shape factor. The assumed boundary layer profile used 
in the present code is a Coles-type profile modified as suggested by Walz [161. 

u - u T [<-3y*-5)e' 3 »* + -^ln(l+y*) + C+ -^w(|-)] (3.25) 

In Eq. (3.25) C is a constant taken as 5.0, II is a parameter which determine the 
strength of the wake component of the velocity profile, and v is Coles wake function 

[16]. 
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Following [4], the radial component of velocity in the fan or turbine strean is then 

given as (3*26) 

v = u a i 

for the fan stream and 


fan 


v = u a 


(3*27) 


core 


for the turbine stream where a is defined as a mean flow angle. 

Finally, the swirl velocity is specified in terms of a mean swirl velocity 
magnitude and by specifying the radial limits over which there is to be swirl 
as an initial condition. The swirl velocity at the initial station (but only 
at the initial station) is assumed constant in the azimuthal direction. This 
condition could be changed at an arbitrary inlet swirl profile, if required, 

but is adopted here for convenience. 

The majority of calculations to be performed with the present 
code are expected to be run without resolving the hub and casing boundary layers. 

For those cases, the initial velocity field in the h-odiate vicinity of solid surfaces 
oust be treated io an approrinate manner. This is accomplished by usrng uall 
functions vhich assume that locally betueen the uall and the first grid point 
auay from the uall the velocity profile is logarithmic soch that 


u . Ur( _L h *L + c) 


(3.28) 


where u x is the friction velocity, < = 0.43 is the von Karman constant, and 


c is taken as 5.0. It follows that 


du 

Am 


U 


r 

rV 


(3.29) 


Equation (3.28) is used to determine a wall slip velocity such that the finite 
differenced form of the velocity gradient one point off the wall is consistent 
with the assumed law of the wall profiles. When wall functions are used to 
represent the primary flow velocity, the secondary flow velocities are still 
required to obey no slip conditions. Thus in the near wall regions the 
secondary flow velocity is assumed to vary from zero to its free stream 
value in a cubic manner with distance from the nearest wall across some 
specified secondary flow boundary layer thickness. Once the velocity and 
temperature distribution are set up, the density is computed to be consistent 
with the temperature and local static pressure as determined from the mean 
static pressure corrected for the normal pressure gradients. 
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lnitia l conditions for the turbulence model arc deter.iued throng), 
specification of a reference length scale, l re£ and free * 
levels. hear the hnh and shrond the length scale is detemn-ed frm. the 

McDonald-Camarata [17] relationship 

<y 


1 - 


(3.30) 


ohere y is the distance free, the surface in question. In th. regm. t 

iohe surface, the distribution of Eq. (3.30) nay be used with , representing 
the distance to the lobe from the grid point in question. H»ever, since 

ss arbitrarily close to a grid point the length scale determined 
“ Z 3-30) will approach 0.0 as y approaches 0.0. This results in^lou 
Reynolds number region in th. initial flou field uhich can create problems 
JL the turbulence model. This problem is circumvented through - of an 
option uhich renders the length scale constant in the region near e , . 

and the distribution of Eg. (3.30) is used only near the hub or casing if 
b darv layers are specified there. With the length scale distribution knoun 
t°b“e initial turbulent viscosity is obtained from the generalited oiling length 

relationship 


rff 1 e fLiif + (1 * (4- 

Pl'P‘ KlT 2 d7) + lh a d,> 'Mz/ 1 hs dir 


(3.31) 


the initial turbulence kinetic energy and dissipation are then detemeined from 
simultaneous solution of Eqs. (2.25) and (2.26). The kinetic energy is 
modified slightly f.llouing McDonald and Ereskovsky 181 to accou 
free stream turbulence, and the turbulent viscosity is then recomputed 

V nnd c to be consistent. 

Under certain conditions, neither the .all length scale option E,. 3. 

nor the make length scale option (constant length scale) , dU 

factory initial values of k, e and ly This may occur ene 

. r n (3 31 } vary greatly between adjacent gri 
rhf> velocity gradients in Eq. vary g y 

points. Under these oitcnmstnnces, it is osefnl to initiate the thence 

quantities using .n option uhich specifies the turbulent kinetic «ergy, 

dissipation and viscosity as essentially constant throoghootthe fl« ■ 

Such on approach has also been used previously hy launder, se. 

Spalding (11! in the computation of free shear flous. 
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If experimental profiles are read la, the specification of initial conditions 
is similar. The velocity components and static temperature are input at each 
grid point. The lobe shape is specified as in the automated starting procedure, 
and a reference length scale and free stream turbulence intensities must be 
provided. Additionally, estimates of the boundary layer thickness on the hub, 
shroud and lobe mast be specified for use in setting up the initial turbulence 
quantities. Initialization of the turbulence quantities is similar to that 
used in the automated procedure. 

Boundary Conditions 

To march the solution downstream, it is necessary to apply boundary 
conditions in each transverse plane. Boundary conditions are needed for the 
energy equation, the primary flow momentum equation, the streamwise vorticity 
equation, the sector and scalar potential equations, and the turbulence model 
equations. The boundary conditions will in general depend on whether swirling 
or nonswirling flow is considered, and whether or not wall boundary layers are 
considered. Referring to Fig. 4, the boundaries at 0^ and 6^ represent 
either symmetry planes or periodic surfaces. Boundaries at the hub (or center- 
line) and shroud surfaces require conditions appropriate for solid surfaces 
(hub or shroud) or a centerline (symmetry). 

Honswirling Flows 

For nonswirling flows boundary conditions along 9 mia and 0^ are 
specified as follows: 

1) aximuthal gradients of U^, E, k and c are set to zero 

2) the vorticity is set to zero, and 

3) the normal velocity component is zero, thus the vector 
potential is zero and the azimuthal gradient of the 
scalar potential is zero. 

Boundary conditions on the hub and shroud surface are specified as follows: 

1) the normal velocity is zero along these surfaces, and the vector 
potential and the normal gradient of the scalar potential are 
thus set to zero, and 

2) these surfaces are taken as adiabatic, and thus the normal 
gradient of enthalpy is set to zero. 

For the majority of calculations it is not expected that the wall boundary 
layer will be resolved. Thus, the appropriate boundary condition for the 
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primary flow velocity, U p , is the gradient condition of Eg. (3.29). For the 
turbulence quantities k and e boundary conditions consistent with the log 
law are 



1/2 (3.32) 

c 

u 


which is used to replace Eq. (3.31) one point off a vail and then 


3k 

3y 


0 


(3.33) 


is specified on the wall. For e the condition 

u 3 

e = _L (3.34) 

icy 

is used. The use of Eqs. (3.32 and 3.34) effectively sets the turbulent viscos- 
ity one point off a wall consistent with the log law. If boundary layers are 
resolved, then Up = k = e = 0 on solid walls is specified. 

With regard to the boundary conditions for the secondary vorticity, 
in Ref. 1 it was argued that a zero vorticity boundary condition would be 
suitable for use with the wall function approach to the primary flow. 

Although this was satisiactory in many cases, in some instances it gave rise 
to large tangential velocity components on the hub or shroud surfaces. As a 
result, under the present modifications of the code the secondary flow vector 
potential and vorticity equations are solved using a coupled solution procedure. 
This procedure allows specification of the correct coupled stream function- 
vorticity boundary condition which ensures that the secondary flow velocity 
components obey no slip on solid walls. 

Swirling Flows 

With the exception of the vector potential, boundary conditions on the 
hub and shroud surfaces remain unchanged for swirling flow. Special treat- 
ment must be given to vector potential, however, because flow may pass in and 

out of the computational region along the radial lines at 9 . and G 

min max 

The amount of flow crossing these radial lines is due to the swirl velocity, 
and a means of determining the growth or decay of the swirl must be accounted 
for. This is accomplished by solving a 9-momentum equation along one of the 
radial coordinate lines prior to solving the secondary flow vorticity-vector 
potential system. The resulting prediction of the 0 velocity component can 
be integrated to obtain the difference in the stream function between hub and 
shroud boundaries due to swirl. 
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(3.33» 



which provides the necessary information for the boundary conditions* In 

Eq, (3*35), w. is subtracted from w, consistent with Eq* (2*3), since the 
9 

velocity predicted from solution of the 0 momentum equation contains both 
rotational and irrotational components* With regard to the boundary conditions 
7 in the azimuthal 0 direction, all flow variables are assumed to be periodic. 

i * 

x 
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Several test calculations were performed to demonstrate the extended 
capability of the modified computer code. The test cases were all performed 
for the geometry shown in Fig. 2 immediately downstream of the lobes with the 
lobe shape shown in Fig. 3. A summary of initial condition for each case is 
given in Table I. As can be seen in Table I these test cases were run with a 
variety of conditions to demonstrate the effect of secondary flows of various 
magnitudes, the effect of swirl, and the effects of the various turbulence 
models. 

Before discussing results from these calculations, it is useful to 
define average values of total pressure loss coefficient, ideal thrust 
coefficient, and Mach number at each axial location, to aid in the evaluation 
of flow properties. A mass-averaged total pressure loss coefficient Cp^ is 


defined as 


r 2(P Tref- P T ) - - 

L />u ' dA 


"PT 


f p\i * dA 


(4.1) 


where ^ is a reference total pressure associated with the turbine strean 
reference conditions, and P re p and u ^ are the reference density and velocity 
associated with the turbine stream. The value of C pT at the starting plane is 
usually not zero due to the presence of shear layers in the initial conditions, 
and for clarity this starting value is subtracted from C pT in the presentation 
of results. The distributions of Cp T presented thus represent losses incurred 
downstream of the starting plane. An ideal thrust coefficient is defined and 
based on the thrust which would be obtained by isentropic expansion from local 
conditions to a predetermined exit pressure. The ratio of the local value of 
this thrust T to the value T^ based on the initial flow field is defined as 
the ideal thrust coefficient, T/T^. The thrust is obtained from the relationship 


T - / r p . M / dA . 


(4.2) 


where is the assumed exit pressure, M g is the Mach number based on 

isentropic expansion from the local pressure to P , and A is the effective 

e e 

exit area, again based on isentropic expansion to P^. 
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Effects of Secondary Flows 

The effects of the magnitude of the initial secondary flows prescribed at the 
lobe exit was investigated by performing two calculations designated as case I and 
case II in Table I, where the initial conditions are summarized. In the case I 
calculation, a radial inward velocity equal to 13% of the local fan stream velocity 
was specified. In the turbine stream, a radial outward velocity of 8% was specified. 
The case II calculation had fan and turbine streams radial velocities specified as 
20% and 25%, respectively. Both calculations used the new optional k-£ turbulence 
model, thus the calculations also serve to demonstrate this capability. 

The relative effects of the different initial secondary flow on mixing are best 
examined by looking at the total temperature contours at the nozzle exit plane. As 
noted in Table I, the temperature ratio between the streams is different for cases I 
and II, however, in the absence of initial secondary flows the thermal mixing between 
the two streams results in a ’tear drop’ total temperature pattern (c.f. Ref. 3). 

The departure from this ’tear drop’ pattern thus indicates the degree to which in- 
itial secondary flows effect the mixing process, even though these two cases have dif- 
ferent initial temperature ratios. The total temperature contours are shown in 
Fig. 5. The case I results show a small but significant departure from the ’tear 
drop’ total temperature pattern observed when no secondary flow is specified at the 
initial plane. However, the results for case II, where even larger secondary flows 
are specified at the initial plane, show a much more dramatic distortion of the total 
temperature contours at the nozzle exit. Here the hot fluid in the upper portion of 
the lobe is effectively pinched off from that in the lower portion of the same lobe, 
and a classical mushroom-shaped pattern (c.f. Refs. 4 and 5), indicating the pres- 
ence of strong counter-rotating vortices is observed. These results clearly dem- 
onstrate the important effects of vorticity or secondary flows in the initial condi 
tions, and indicate the importance of lobe generated secondary flow on the tempera- 
ture profile development. 

Effects of Turbulence Models 

The influence of the various turbulence models, k-e, k-£ or the wake model 
were investigated by rerunning case II with the k-c and wake turbulence models as 
cases III and IV, respectively. With the exception of the turbulence model used, 
cases II, HI and IV are identical as can be observed in Table I. The results of 
these calculations are summarized in Table II where area-averaged values of total 
pressure ratio, and total temperature ratio are presented along with the previously 
defined total pressure loss coefficient, and thrust ratio. Examination of the infor- 
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nation presented in Table II indicates that the computed results for this lobe mixer 
flow show only very limited variation with the turbulence model. The results indi- 
cated that the k-e model gives the higher losses whereas the wake model predicts the 
lowest. In view of this limited dependence on turbulence model, it would appear that 
turbulence modeling in lobe mixers is not an area of primary importance. This is 
born out further by the comparisons presented in Fig. 6 where the total temperature 
signatures resulting from each calculation are shown. This comparison shows that 
only small scale details are altered by the different turbulence models while the 
large scale features remain unaffected. Additional evidence that turbulence modeling 
may not be of primary importance in lobe mixer flows is provided in Fig. 7 where 
the secondary flow vorticity contours at the nozzle exit are compared for the 
three turbulence models. Here also, it is observed that the predicted vorticity 
contours for all three turbulence models are nearly identical. The pattern of 
two large counter-rotating vorticities is observed for all three turbulence 
models, and the strength and size of these vorticities remain uneffected by 
the turbulence model used. Some small differences in the contours are observed, 
however, particularly when using the k-i. turbulence model. These small scale 
differences appear as slight distortions of the vortex contours, but as indicated 
in Table II, this effect does not apparently alter the mixing in any significant 

manner • 

Effect of Swirl 

Perhaps the most significant modification to the computer code performed 
under the present effort is the extension of the capability of the code to 
swirling flow. To demonstrate this capability several computations were carried 
out for swirling flows. The first of these calculations, case V, is identical to 
case III discussed previously with the additional specification of an initial 
swirl component of velocity, equal to 10% of the hot stream velocity, across the 
entire mixing duct (hub to shroud). The pertinent integrated quantities for 
these two calculations are presented in Table III. 

As would be expected, case V with swirl exhibits higher losses reflected by 
the area averaged total pressure ratio and the total pressure loss coefficient, 
and a lower thrust ratio. The contour plots of the exit plane total temperature 
distribution for these two cases is shown in Fig. 8. This comparison demonstrates 
that the large scale features of the signature remain almost unaltered by swirl 
and are simply convected in the azimuthal direction. Some distortion of the 
pattern is observed, particularly as the center of the mixing duct is approached; 
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however, these effects do not appear to enhance the mixing of the two streams 
significantly. A comparison of the vorticity contours at the nozzle exit plane, 
shown in Fig. 9, indicates a more substantial distortion of the counter rotating 
vortex pair by swirl. This suggests that in longer mixing ducts, where swirl 
effect may develop further, more significant effects on mixing may be observed. 
However, in view of the higher losses incurred with swirl, and the absence of any 
clear-cut benefits resulting from swirl, it would appear that swirl in lobe mixers 
may not be beneficial from a performance standpoint. 

A final calculation with swirl, case VI, was made using the wake turbulence 
model, and with 10% swirl velocity confined to the outer half of the mixing duct. 
This calculation was made to determine if residual swirl in the fan stream would 
alter the mixin g significantly. The exit plane total temperature contours for 
this calculation are presented in Fig. 10. Again, no significant effect on 
mixing is observed. The large scale characteristics of the total temperature 
signature remain intact with some distortion observed in the outer region of the 
flow. 

Although swirl did not appear to have a significant effect on the mixing 
characteristics of the lobe mixer/exhaust nozzle used in the present calculation, 
it should be noted that the exhaust nozzle was short. If a longer exhaust 
nozzle were used, the effects of swirl would have more time to develop, and 
significant effects on mixing could result. However, as noted earlier, swirl 
does appear to give rise to higher losses, and if performance is the objective, 
swirl does not appear to be beneficial. 
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SUMMARY AND CONCLUSIONS 


An analysis which was previously developed for detailed prediction and 
computation of three-dimensional subsonic turbulent flows in turbofan lobe 
mixers has been extended and modified. An improved starting procedure which 
allows specification of non-zero starting vorticity using a parametric repre- 
sentation of radial flow velocities at the initial plane has been developed. 

The method of scaling an incompressible potential flow to obtain the compres- 
sible flow imposed pressure gradients has been improved. An optional k-J. 
turbulence model has been incorporated within the code, and the original k-e 
turbulence model has been modified to provide improved reliability. Finally, 
a swirling flow capability has been added to the code. 

Calculations have been performed for a number of test cases to verify 
the modifications to the code, and to demonstrate the importance of the effects 
of secondary flows, turbulence model, and swirl on mixing. The results of these 
calculations indicate that secondary flows generated within the lobes play the 
most dominant role in determining the degree of mixing. The effects of the 
turbulence model appears to play a very limited role, altering only small scale 
details of the predicted flow field. Finally, calculations where a swirl component 
of velocity, equal to 10% of the turbine stream velocity, was specified across the 
entire mixing duct indicated that swirl may not alter the degree of mixing sub- 
stantially. In this case, the major effect was to simply rotate the total 
temperature pattern observed without swirl with little additional distortion. 
However, specifying swirl over a limited portion of the duct radial extent may 
have a more pronounced effect. More detailed investigations and experimental 
confirmation will be required to fully assess these effects. 
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CASE tf 


lp/U T 

V t t 

V F 

.633 

.4 

-.13 

.86 

.74 

-.20 

.86 

.74 

-.20 

.86 

.74 

-.20 

.86 

.74 

-.20 

.86 

.74 

-.20 


Up “ Fan stream velocity 
F 

U T ■ Turbine stream velocity 

T- ■ Fun stream temperaturo 
F 

T » Turbine stream temperature 
T 

Vp ■ Fan stream radial velocity 
F 

V * Turbine stream radial velocity 
T 


i 


TABLE 1 


w s i. 


k TURB MOD 


0 .006 .04 


0 .006 .04 


0 .006 .04 


0 .006 .04 


,25 .1 full .006 .04 


outer 

.25 .1 half .006 .04 


VI ■ Swirl velocity 
s 

jt « Turbulence length scale 
H Turbulence kinetic energy 


Re « Reynolds number 


TABLE II 



TURBULENCE 

TOTAL PRESSURE 

TOTAL TEMP. 

C PT THRUST 

CASE 

MODEL 

RATIO 

RATIO 

LOSS = (+) RATIO 

II 

k -l 

.9892 

.7866 

.0956 

.9930 

III 

k-e 

.9890 

.7863 

.0975 

.9926 

IV 

WAKE 

.9896 

.7866 

.0923 

.9933 
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TABLE III 




CASE 

TURBULENCE 

MODEL 

TOTAL PRESSURE TOTAL TEMP. 
RATIO RATIO 

C PT 

LOSS = (+) 

THRUST 

RATIO 


Ill k-e .9890 .7863 .0975 .9926 


.973 


.7884 


1132 


9883 
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Figure 1. - Schematic of coordinate system for 
three-dimensional flow problem. 
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Figure 2. - Typical lobe mixer geometry and potential flow coordinate system. 
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ect of initial secondary flows on 
.1 contours. 


exit plane 
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Case III 
No Swirl 


Case V 
10% Swirl 



Fig. 8 - Effect of swirl on exit plane total temperature contours. 




APPENDIX - USER'S GUIDE 

Introduction 
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The Mixer computer code may be run using initial conditions obtains 
fro. experimental data and read in on a pointuise basis, or an initial no 
field my be generated by the code itself open p.ra.etric specification o 
the proper inpot variables. In either case, the originally spec, led fov 
fields ate processed by the code to assore that the initial flow field is 
const tent with the finite-difference approximations to the governing egoa.ion. 

-• — ■ u is a,s :.r::; 

to activate a swirling flow option of the .ode. Worn this option is util 

, „ , , . ..i-indic ;n»! t initial conditions, with the 

j»„v flow is assumed to be period it, 

exception of the swirl velocity, arc assumed to be symmetrical about the lob • 

center I ine. Calculations perforned with swirl must, therefore, const era u 

il fin* verse planes, whereas only a half lobe nee 
lobe within the computational tran.tver.se. j. 

considered for nonswirling flows. 

an ot fort to keep the initial coodifion specif icat, on st.tlar for both 
... lrHll8 and nonnwirl ing flow calculations. the initial flow field for swirling 
other tn, co, a, lone performed os, eg periodic boundary coadtttnos, as ... 
half the i.,be and then reflected, a a result, the user may use the same 

: P „, rtm stream fc both periodic and non, , t iodic calculations by ~ 

„f .-rid points, the asinuthal cor.potat tonal limits, and calling for 
options. Specif leal lor. of the 1 be shape need not be altered; only 

ba , • ' a lobe must be specified via the input in either case. 

Th ,;,, r , h cuid read over all input descriptions before attempting to set up 

a case. 



description of input 


Input to the MIXER code is accomplished by a combination of card input 
and information which is stored on auxiliary files. Metric information 
generated by the ADD code [13] and restart data sets are stored on separate 
auxiliary files. Program control is performed through card input, which is 

divided into three categories: 

1) Plot file information for NASA— Lewis plotting package, 

2) NAMELIST input, and 

3) Experimental flow field input. 

Format of the required input is as follows: 


Plot File Input 


CARD # 

COLUMNS 

VARIABLE 

Card 1 

Col. 1-32 

TITLE (I) 

Card 2 

Col. 1-2 

ISYM 


DESCRIPTION 

Plot title, format 5A6.A2. (May be 
left blank) . 

o 

Twice the number of lobes in a 360 
cross section of the mixer geometry. 

Indicates the type of coordinate system 
for plot routines. Format F10.0, 

Use SYSTEM=2. 


Col. 3-12 


SYSTEM 



NAMELIST INI 
Namelist 

$IN1 

RG 

CP 

tzero 

13ZERO 

PZEROl 

TZERO 

VISC1 

$END 

NAMELIST IN2 

The desired 
namelist $IN2. 

Namelist 

$IN2 

IAXI 


XENTR 


NAMELIST INPUT 


Description 

Cas constant, ft lb £ /(ulu 8 °E), default 1716.3 (air). 

Constant pressure specific beat, ft lb £ /(sl»t °M. 
default is 6012.384 (air). 

Reference temperature, assumed to be nominal primary 
stream temperature in degrees Rankrne. 

Reference velocity, assumed to be nominal .primary 
stream velocity, ft/sec. 

Reference static pressure, f^sumed default is 2116-8. 
pressure at initial plane, lb f /ft , dera 

Reference length, ft. 

Reference viscosity, slugs/ft-sec. 


computational options 


and mesh information are specified in 


Description 


=1, cylindrical coordinates. 

=2, rotated orthogonal coordinates generated from 
ADD code geometry file. Default. 


If IAXI is input = 2 
XENTR, ILAP, LX, LY, 


(or defaulted) the variables, 
and IPCOR must also be input. 


Dinensiouless dlutancc, fru. lultlul - 

line to initial plane of lobe mixer calculation. 

Used only for IAXI=2. Determined from ADD code 
parameters as XENTR=(JSTEP XENTR -D*DS where JSTEP XENTR 

is the number of the ADD code potential surface 
representing the MIXER code initial plane and DS is the 
streamwise delta used in the ADD code. See Fig. Bl. 
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Mamelist 

ILAP 


LX 

LY 

IP COR 
ICOMP 

IDATA 

ILAM 

IDIS 

LOWRE 


Tiggcriptioo 
= 1 for IAXI=2 

= 0 if IAXI=1 If IAH is set to 1, 

ttap is defaulted - 1- K 

Ut is the »«■*«' ° £ v °^,st he inpet i £ Ih*M- 
code geoaecry ^ the c«ie 

« is the nuaber of^^ ^ lnfat u TUX-1- 
geometry • „ UMtlfy the potential 

IPCOR is used with ILA ^. ensiot iless velocity is 1-0 
surface whose average diuensx iPCOR 

in the axisymmetr -rSmu-TSTEP in 8 eneral * 

should be input as IPCOR-JSTEP^^ 

= -1, stop after input dump- 
- n normal run (default). 

I ,’ Ip after printing initial Plane geometry. 

I 2 ’, sto p after printing initial flow field. 

= o, automated start (default). ‘ 

, „*** and temperature iiexa 

= 2 read initial velocity and temp 

d.u at each grid po»t. 

ILAM=1, Laminar flow. 

ILAM=0» k-e turbulence model. 

ILAM--1* Wake turbulence mode . 

ILAM=— 3, k-t turbule “ ce rb “^ence kinetic energy, at 

* SbS s tr eanTvalu e throughout the initial 

=0, initial dissipation, 

turbulence productio q 

TTVTS=1 is recommended. • 

=0, use high Reynolds number form of tur u en 
model equations 

-1. use low Reynolds number «££££ mo^ ^ 
equations. Requires IDISJ, ^ Also requires 

specification of wal1 1 ° i(m near wall regions 
a high degree 0 f resA ± appropriate 

using generally recommended 

EPS(I) • Th is option is nor g 
for mixer flows. 
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Namelist 


Description 


YS(I, J) 


NE(I) 


IGRID(I) 


EPS (I) 


XCTR(I) 


NS 

AP 


YS(I, J) sets the limits of the transverse computational 
planes. 

YS(1,1) = r lir - n (dimensionless, r/YZERO) 

YS(2,1) = r Tnax (dimensionless, r/YZERO) 

YS(1,2) = 0^ (radians) 

YS(2,2) = (radians) 

For mixer geometries generated by the ADD code YS(1,1) must 
be the dimensionless plug radius at the first ADD code sta- 
tion RMIN, and Y(2,l) must be 1.0. See Fig. R-l. For non- 
swirling or non-periodic flows, YS(1,2) must correspond to 
the 0 value at the maximum lobe penetration and YS(2,2) to 
the 0 value midway between the lobes. See Fig. B-2. For 
periodic or swirling flows, YS(2,1) and YS(2,2) correspond 
to sequential values of 0 midway between the lobes. See 
Fig. B— 3. 


NE(1) = number of grid points in the radial direction. 
NE(2) = number of azimuthal grid points. For periodic 
calculation NE(2) must be odd . 

Maximum values for NE(I) are 50. For periodic calculation 
NE(2)max = 49. 


IGRID=0, no grid stretching (default). 

=1, Roberts stretching at r min or 0^. 

=2, Roberts stretching at both r ^ and r MX 

or 6 . and 0 • 

min max 

=4, hyperbolic sine stretch about XCTR(I). 

1=1 for radial direction, 1=2 for aximuthal direction. 


Used for IGRID(I)/0 to control the grid stretch near bound- 
aries. For Roberts stretching 0.0<EPS< f.0. For hyper- 
bolic sine 1.0<EPS<5.0. 


Used only with IGRID(I)=4 to specify the r(I=>l) or 
0(1=2) location about which grid points are clustered. 
YS (1,1) < XCTR(I) < YS(2,I). 


Number of axial stations to be computed. 

NS=1 is initial plane. Maximum value is 50. 

If AP is input .GT.O a geometric streanwise grid is 
set up with Ax(J)/Ax(J-l) = AP. 


If AP is input .LT.O the user must specify the axial 
mesh. AP has a default of 1.05 but the AP.LT.O option 
is recommended. 
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X(JX) 


IPERIO 

ISWIRL 

Send 

NAMELIST IN3 
Restart 
Namelist 

SIN3 

IRSTIN 

IRSTOT 

JRSTIN ' 
JRSTOT 


X is an array which contains NS values of the stream- 
wise mesh coordinate. 

If AP is input .CT.O only X(l) and X(2) need be input. 

If AP is input .LT.O NS values of X must be input. 

If IAXI=1* X is a physical but dimensionless distance. 

If IAXI=2 and ADD code geometry is used, X represents 
“.^iLsionlnns co.nutntfonnl *«»»“ from the 
initial plane. 

The maximum allowable value of X is given as 

(JSTEPhax- jstep xentr)*ds/uscale 

L/, lV-YSa 1))- The maximum value of X should De 
Shiw less than this valne to avoid Mruhnng ont of 
the ADD code geoaetry which will result an an error 
termination- See Fig. Bl. 

=0, normal calculation 

=1 periodic option activated. Input must be set up 
fo^ periodic flow field. Check NE(2) specification. 

=0, nonswirling flow. 

.1, swirling flow. Requires IPEHCl. Also requires 
values for RSU1, RSU2, VSURL, and DE1.SUI. 


and 


plot file information is specified in namelist $IN3. 


Description 


The number of the axial 
restart. No restart if 


station to be read in for a 
IRSTIN=0 (default). 


Increment for saving restart information. Restart ” he 
be written every IRSTOT steps from starting value of JX, th 
axial station counter. Default assumes no restarts writte . 






File number from which restart information is to be 
read. Default file number is 11. 

File number onto which restart information is to be 
vritten. Default file number is 11. 
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Namelist 

NFILE 

NSAVED 

1PL0T 

$END 

NAMELIST IN4 

Velocity and 
namelist $IN4, as 
initial secondary 

Namelist 

$IN4 

USTRM2 

TSTEM2 

TUSTRl 

TUSTR2 

ALEN 

amexit 


inscription 

Sequence number in file 
information. 


JRSTIN cf the desired restart 


Xbe number of restart data blocks saved o» file 
JRSTOT . 


If IPI.0T21 a plot file oust be 

stream and information " l IPLOT"0 gives no 

fnr subsequent use in plotting- 

Plot file! Default is IPLOT=0. A plot foie will be 

* nn 


temperature ratios between the two streams are specified in 
are turbulence levels and information used to construct the 

flows. 

Description 


Ratio of secondary stream (fan) velocity to prinary 

stream (tnrbine) velocity. For matehedj 

Mach numbers and total pressures USTRM2- ,/TSTRMZ. 

Default value is USTRM2=. 707 . 

Ratio of secondary stream temperature to primary stream 
temperature. Default is TSTRM2-.5. 

Primary stream turbulence intensity- 

X = u f ^/n? where u=UZER0- 
x u± 1 i 

Secondary stream turbulence intensity- 
X = u* ^/ u 2 U 2 = ^Z£.R 0 xUSTEM.£ - 

If 2 TUSTR2 is input 0.0 the free stream turbulence 
in both streams (k=3T u 2 /2) is assumed equal to 3T U1 /Z. 

If IDIS-1 is input (as recommended) only TUSTRl is used. 

Dimensionless free stream turbulent length scale, 
Iref/YZERO. To be estimated by the user based on 
ge^meSc constraints of the mixing duct and boundary 
lavers. ALEN should be on the order of .1 of 
estimated free shear layer or boundary layer thickness . 

AMEXIT is an assumed exit plane Mach number which 
S to specify the iseuttopic exit plon. press«e for 
the thrust calculation. Default value xs AMEXIT 0.9. 
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Namelist 


Description 


DELMIX 

DELHUB 

delsh? 

RHUB 

RTIP 

DELSEC 

VFAN 

VCORE 

VSWRL 
RSW1, RSW2 

DELSWL 


DELMIX is input with IDATA .GT. 0 as a dimensionless 
boundary layer thickness presumed to be representative 
of the boundary layers on the lobe surfaces. It us 
in computing the initial turbulence length scale. 

DELHUB is input with IDATA .GT. 0 as a dimensionless 
layer thickness, 6/YZERO, respresentative of the bub 
layer. It is used in computing the initial turbulenc eig 
scale. See IHUB in $IN5. 

Similar to DELHUB but pertaining to the shroud boundary 
layer. See ISHR in $IN5. 

RHUB is a dinensionaless secondary flow hub bou "^ y 
“^thickness used to eonsttoct initiol secondary 
flow profiles required for IDATA- 0 or 2. 

Similar to RHUB but pertaining to the shroud secondary 
flow boundary layer. 

A vaSe on the order of (RTIP+RHUB) /2 should be used. 

Used with IDATA=0. VFAN is the ratio of *** 
velocity to the local streamwise velocity (positive 
outward) in the fan stream. 

Used with IDATA=0. VCORE is the ratio of the local 
radial velocity to the local streamwise velocity 
(positive outward) in the turbine stream. 

VSWRL is used only with ISWIRL=1 to specify the 
dimensionless swirl velocity magnitude. 
n<;ui and RSVP are the radial limits over which swirl 
ITLZ 5 be specified. They .use be input such 
that YS(1,1) < RSW1 < RSW2<-YS(2,1). 

OEbS f is 

value. V6 DELSVL must be less than (RSW2 - RSWl)/2. 


$END 

“the shape and primary flow boundary layer parameters are 
specified in namelist $«5. Specification of this information should 
performed using Fig. B2 or B3 as a guide. 


Namelist 

$IN5 

ILOBE 


Description 

ILOBE is a flag used in the automated starting routine. 

ILOBE=0 indicates that no boundary layers are to be 
set up on the lobe surfaces. 

ILOBE 3 ^ indicates that lobe boundary layers are to be 
set up. 


Namelist 


Description 


IHUB 


XSHR 

bliiub(i,J) 


NLOBE 

YLOBE(I.J) 


present in the u input .CI.O end IHUB is 

starting options. If If the IDATA option 

input .GT.O DELHUB must be i"^* are des ired, IHUB 
is not used, and hub bound ^ pcints inp ut to 

rscr'ibe'X'hub bounder, leyet in the err,, BIHUB. 

ISHE is similar to IHUB estept it pertains to the shroud 

boundary layer. 

If hub and/or sh roud boundary layers are ^construe 

ted using the automated s ^ r ^ kness and shape factor 
boundary layer disp n ^iens of 6 along the hub and/or 
must be input as u^ accomplished using the array 

shroud surfaces. This is au 

BLHUB. See Fig. B2. 

i no nf G (in radians) where 
BLHUB(I.l) contains the ™£s are specified. The 9 
hub boundary layer pararu s ffiUSt incr ease monotonically. 

values are arbitr ry IP ERIO=Q option the first 

For calculation using t ^ the ^ ast value 

vlloe must bn slightly f 3 ” 8 "“ ,“ t IPEEIO-1, the 

firsX XfeXsfbfsUg^tly 

Sf^e^f «‘v*^“Xeci£ieH must be n,ual tn the 
input value of IHUB. 

. fl * /V7FR0 the dimensionless displacement 

“SSss’rf the hd, boundary layer at each speerfred 
0 location. 

ir or pach 9 location. 

BEHUB(I,3)=H, the shape facto 

. , oniri valent of BLHUB(I,1)- 

BLHUB(I,4) ^ the shroud ^ BUIUB(I,4). 

ISHR values of 0 must d 

BLHUB(I,5) and BIHUBU.M ate similar to 
and BIHUB (I, 3) except they pertain 
boundary layer parameters. 

NT QBE is the number of data points input to describe 
Te lobe shape and lobe boundary layers. 

YLOBE is an array ^ hich ^boinda^iayers^Sch may 
specify the lobe shape an ^ NLOBE values of 

be on each aide of the lobe surfac^^^^ ^ 2 

each parameter in specification of the lobe 

XpfXXX "Xiob?ctnsses the ^ bnnndsty. 
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List Descri£tion , ess radial coordinates, 

YLOBEd, D contains ^^ “^st point must be the radial 
E(I>J) r/YZERO, of the lo J e - ? e crosse s the 6^ boundary. The 

coordinate where the J ^ d f al coordinate value where the 

££s --- «*. ♦ *->«• 

See Fig* R va i ue s 

EtOBE(I.» c «nc. ? s ^values correspon^ .^lightly 

r°izT". — 
H 2 SSS i e e“ S5£ ♦ 

«“St te slighC 7 dimensionless dis- 

turbine stream si YLOBEd.D and YLOBEd, >* 

locations specill perta ins to 

rhe sa me as YLOBEd, 3) a* 061 * 

the^f an* stream side of the lobe. ^ ^ ^ factors> H, 

a t^ each* coordinate 3 point^ f or* th^ tur b in e stream side of the 

i„. ~ ••••• 

for Che Ee» stt “” ' , £irf 1£ thc 

EL 0 BE(I. 3 ) through YLOBE(E»B) need ItW » and 
IDBTA option IS of the lohe, -* 

YLO e B cifi^ fot all starting options. 

SP . „ f ac tor in BUTOB and EhOBE 

If IWMCE-1 the slots th e uafce function 

IUAKE need “LSdaTlayaP P- £ile 15 aSS ” ed ' 

SEtm taeeous to consider specification of 

hhen setting up . loBe shape it is a van th ese conditions 

bl „BB and ELOBE for a ‘ „ g . B2 . I£ e„ ln is specified as rem. 

the computation segment appears ^ i§ pI „perl, specified. 

and i. specified . s Posltive, changing the E«»E 

then a periodic ealcula . .6 , adjusting the mui4.er of 

or BIHDB arrays simply by setting min »a* 
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V 


e grid points, and by specifying IPERIO-1. In general, if the arrays YLOBE 

and BLHUB are properly specified for a nonperiodic calculation of limits 0 min 

to a the same lobe shape specification will be correct for periodic calcula- 
to °max * 

tions in which 9 may remains unchanged, and 6 min is reset as 6 mi n = ^®min -0 max‘ 
NAMELIST IN 6 

Boundary conditions are specified in namelist IN6. For periodic flows 0 
direction boundary conditions are ignored. The boundary conditions are set to 
defaults appropriate for mixer calculations performed with wall functions. 

Only XCBC must be specified. 


Namelist 

$IN6 

NBCON(I,J,K) 


NBC0N(1,J,1) 


NBC0N(2,J,1) 


NBC0N(4,J,1) 


Description 


NBCON is a three-dimensional array which allows the user 
to select the desired boundary conditions for the primary 
flow momentum equation, the energy equation, the secondary 
flow stream function-vorticity system, and the turbulence 
model equations. The subscript "J" determines the 
boundary surface. J=1 for the surface at for 

R max» J=3 for 6 min» for 0max ’ 

NBC0N(1,J,1) is used to specify boundary condition for the 
primary flow momentum equation. 

= 0 - no slip 

= 1 - zero normal gradient 

= 2 - zero second derivative 

= 4 - wall function 

At the hub (J=l) and shroud (J=2) either no slip or wall 
functions should be specified. (Wall functions are preferred, 
and required when L0WRE=0 in $IN2) . 

At 0 . (J=3) and e„_„ (J=4) symmetry (zero gradient) 

min max 

conditions are required when IPERIO=0. 

NBCON (2, J,l) is used for the energy equation. 

= 1 - zero normal gradient (adiabatic) 

= 2 - zero second derivative. 

Adiabatic conditions should be specified at all surfaces 
(J~l to 4) . 

NBC0N(4,J,1) is used to specify boundary conditions for 
the vorticity-vector potential system, and results in 
either no slip or slip conditions applied to the local 
tangential component of the secondary flow velocity. 

* 0 - no slip 

« 1 - slip with zero vorticity 
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kJ 


NBC0N(1,J,2) 


NBCON(2,J,2) 


At r (J=l) and r (>2) no slip should be specified. 

a m in / an( j e max (J=4) symmetry implies slip with 

“ o^rSt, specif- The po slip ondici.. 

should be used at a centerline. 

tIBCON ( 1 J,2) Is used to specify boundary conditions for the 
3S kinetic energy equation uhich is solved if 

au O 


= 0 — zero kinetic energy 
— 1 — zero normal gradient 


= 2 — zero second derivative 

... _ / t=i'j and r (J=2) the turbulence kinetic energy 

shoul^be set to ze%if L0WRE=1, or the normal gradient 

should be set to zero if LOWRE-O. The latter is consistent 

with the wall functions used for the primary flow momentum, 
with the wan and 0 m _„ (J=4) zero normal 


and is recommended. At 0min ^ ® 
gradient should be specified. 


D max 


NBC0N(2,J,2) is used to specify boundary conditions for 
hKo t-iirhulence dissipation equation (ILAM=0 only). 


= 0 - zero dissipation 


= 1 - zero normal gradient 

_ 2 - zero second derivative 

= 4 - wall function type boundary condition 

At r (J=l) and r (J=2) the dissipation should be 
set to n zero if LOWRE^f, or the wall function should be 
used if LOWRE=0. The later is recommended. At 0 min an 

B zero gradients should be applied, 
max 


XCBC 


CCBC is used in MIXER calculations to change boundary 
ronditions as the flow moves downstream off the hub, 
/here wall function boundary conditions are applied, 
to a center line where symmetry conditions are applied, 
rhe user must specify the value of XCBC by examination 
if the geometry from the ADD code. XCBC is com puted 
is XCBC=(JSTEP XCBC - JSTEPxeotr)*^ where JS^XCBC is 
the number of the ADD code potential lin e w here the 
plug degenerates to a center line and JSTEPxqjjr is the 
number of the potential line representing the Initial 
MIXER calculation plane. See Fig. Bl. The default 
value causes no changes. 


56 




NBCONC(I 


$END 

NAMELIST 

$IN7 

IPRN(I) 


ICON 


f j t K) NBCONC has the same function as NBCON but is used to 

specify which boundary conditions are to be changed 
for X>XCBC and what they are to be. Default value 
will cause no change in boundary conditions regardless 
of the value of XCBC. For a calculation with hub 
boundary layer, NBCONC is used to change from wall 
functions to symmetry conditions as the hub degenerates 
to a center line. 

IN7 


1PRN is an array which selects variables to be printed 
at each axial station. If IPRN(I)=1 variable "I" will 
be printed. If IPRN(I)=0 variable "I" will not be 
printed. The variables are numbered as follows: 

1= Variable 

1 u 

2 v 

3 w 

4 p 

5 E 

6 

7 k 

8 c 

9 S 

10 * 

11 M 

12 static temperature ratio T/T r 

13 total temperature T°/T° 

14 pressure coefficient 

15 static pressure ratio P/P r 

16 total pressure ratio P°/P° 

17 swirl angle, degrees 

18 u’v 1 

19 tT'w 7 

20 v 1 w * 

ICON controls printer contour plots of the printed 
flow field variable. 


IC0N=1 produces printer contour plots. 
IC0N=0, no contour plots. 


IPRINT 


$END 


IPRINT must be input 2l. Printout will be made 
every IPRINT' th station. 


Experimental Flow Field Input 

Th. following cards are needed only 11 starting profiles are to be read 
• nr iata-?^ Each card contains velocity components u, v, w m feet 

in (IDATA-2) . Each car in degrees Rankine at each computational 

per second and static temper^ ^ thus . total o£ NE(1)*NE<2) 

:::: — - - *• - — **“ p ° int 
numbering shown in Fig- B4. 


VElULALJf 

direction, ft/sec. 
velocity component in tl 
ft/sec. 

velocity component in tV 
ft/sec. 

o 

T - static temperature, R- 


CARDS 8 $ENEBT 

to ? 

Col. 1-10 

u - 

11-20 

v - 

21-30 

w - 

31-40 

T - 
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Radius 



Axis 


Figure Bl. - Typical coordinate system and definition of input parameters. 
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Casing boundary layer parameters 

specified along this surface as 

a function of 0 Boundary 2 



0 


Figure B2. - Specification of the lobe shape and hub and casing 

boundary layers. 
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* . + '* 
m i n max 




Figure B3- - Specification of lobe shape for 
periodic flow domain. 


61 









